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Abstract. We examine the current status of the spectral analysis of 
SN1987A during its early stages. Issues of the shock breakout and UV 
flash, the density and chemical structure, and masses of different layers 
are discussed. A decade later, several aspects need a fresh look and inter- 
pretation. We summarize what questions have been answered and where 
some results disagree. Unresolved problems such as the excess of s-process 
elements and influence of asphericity are addressed. Finally, SN1987A is 
considered as a test case for using Type II as distance indicators via the 
Baade-Wesselink method. 



1. Introduction 

The emitted light of Type II supernovae allows study of the explosion mech- 
anism, testing of hydrodynamic models and probing of the final stages of the 
stellar evolution. The observed spectrum gives direct information on the physical 
and chemical conditions of the supernova photosphere at a given time. Deeper 
layers are observable at later times. Spectral analysis provides an effective tool 
to scan through the atmosphere, to reveal the density, chemical, and velocity 
structure, to probe explosion models, to reveal mixing processes and departures 
from sphericity, and to determine the distance. At any one time, the photosphere 
spans only about 10 to 1O _2 M0 and, consequently, good time coverage is crit- 
ical to derive integrated quantities such as the total mass. The following analysis 
is mainly based on our own work but, for comparison, we will refer to results of 
other groups as appropriate (Table 1). Specific references are given in the text. 

Table 1. Some Work on Atmosphere Models of SN 1987A 



First Authors 


Phase 


Topic 


Ensman / Hauschildt 


0-4 d 


shock breakout, early spectra + colors 


Hoflich 


0-220d 


shock breakout, abundances, structure, 
mixing, asphericities, CO, distance 


Lucy/Mazzali 


2-100d 


structure + abundances 


Jeffery 


2-30d 


asphericity 


Eastman 


2-10d 


structure/distance by Baade-Wesselink 


Schwarz 


150-500d 


gamma excitation 


Li/McCray 


200+d 


gamma excitation, mixing 


Wagoner / Chilukuri 


2-50 d 


distance by Baade-Wesselink 
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Figure 1. Observed 
my (with measurement er- 
ror bars) of SN1987A in 
comparison with that pre- 
dicted by a model with 
1.25 x 10 51 erg (Shigeyama 
Sz Nomoto, 1990) using 
NLTE and extended grey 
atmospheres (EGA). The 
distance is taken to be 48 
kpc. E B _ V is set to 0.15 m . 
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2. Model Assumptions and Free Parameters 

Most of the calculations summarized here were performed prior to 1989 and, con- 
sequently, the description of our Nlte code for Extended ATmospheres (NEAT) 
below corresponds to this era (Hdflich 88ab, 90, 91a; hereafter H stands for 
HMich). It has since been updated. The system of equations and their solution 
is briefly described as follows, a) The radiation transport equation is solved in 
the comoving frame and includes relativistic velocity fields and line blanketing 
(Mihalas et al., 1975, 76ab). Time dependent terms can be neglected because 
the radiative time scales arc short compared to those of physical quantities at the 
photosphere, b) The statistical equations are solved by a perturbation method 
(ALI). For the early light curve (LC), time dependent statistical equations are 
solved (H91a). c) The energy equation. During later phases, we assume homol- 
ogous expansion and radiative equilibrium to determine the thermal structure. 
The bolometric and 7-ray luminosities are taken from the observations. The 
7-ray energy deposition is treated following Colgate et al. (1980). During the 
shock breakout, the temperature structure is taken from the hydro because it is 
governed by adiabatic expansion and the corresponding recession in mass of the 
photosphere (see also Hauschildt & Ensman 1994, hereafter HE94). Detailed 
atomic models are used for up to the three most abundant ionization stages of 
several elements (H, He, C, N, O, Na, Mg, K, Ca). A large number of LTE-lines 
is included in a pure scattering approximation. 



3. The Shock Breakout 

SN1987A is the first Type II supernova to have been observed shortly after the 
initial event (e.g. McNaught 1987). Measurements of the envelope mass and 
explosion energy based on the early LC are consistent with the interpretation 
of the late LC (Arnett, 1988; Nomoto et al., 1988; Woosley, 1988); however, 
flux limited diffusion calculations gave discrepancies between observations and 
models of up to 2 magnitudes at the earliest phases. 

To address this problem, H91a coupled the hydrodynamical calculations of 
Shigeyama & Nomoto (1990) with a detailed atmosphere. The progenitor has 
a total mass of 16 M , 10.4 M of H, a radius R p h = 3 x 10 12 cm. Explosions 
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have been investigated with E^ in of 1.0 x 10 51 and 1.25 x 10 51 erg which give a 
reasonable reproduction of the late light curve. Both NLTE and scattering effects 
are important (Fig. 1, Hoflich et al., 1986, H90, H91a, Pizzochero, 1990). The 
NLTE fluxes differ from those resulting from a flux limited diffusion approach or 
the solution for an extended gray atmosphere (EGA), and the slope of the LC 
changes. In particular, Ly shows a strong local minimum shortly after shock 
break out. The general tendencies, and the drop in Ly, can be understood as a 
consequence of the rapid drop in p and T and the resulting size of NLTE effects 
and of the scattering optical depth (H91a). 

Table 2 gives the results of NLTE atmosphere models during the shock 
breakout phase; the electron temperature T, the radius of the photosphere R, 
the particle density N (in cgs) at the Thomson optical depth of 2/3, My, U-B, 
B-V and the bolometric luminosity L^oi based on the hydrodynamic calculation 
of Shigeyama & Nomoto (1990). The dynamical model had logarithmic zoning 
that resolved the photosphere at the epoch of the Jones' limit. The functions 
Cy and Cm give the ratio of the visual and bolometric fluxes, respectively, 
calculated by NLTE models to those of the extended grey atmosphere (EGA). 
Observed colors are in good agreement (t=4.4 x 10 4 s B-V=-0.2 ... 0.0; at 
9.83 x 10 4 s U-B=-0.836, B-V=0.085, Moreno & Walker 1987). 



Table 2. Results of NLTE calculations of Shock Breakout 



time [s] 


T[K] 


R 


N 


My 


U-B 


B-V 


L(, i 


Cv 




6175. 


788000. 


3.36E12 


6.3E14 


-11.75 


-1.23 


-0.30 


4.6E44 


0.18 


0.39 


6400. 


399450. 


4.20E12 


7.8E13 


-10.99 


-1.28 


-0.32 


8.7E43 


0.13 


0.32 


7000. 


282000. 


6.28E12 


3.2E13 


-11.19 


-1.33 


-0.33 


5.0E43 


0.12 


0.28 


8000. 


102000. 


9.56E12 


1.1E13 


-11.46 


-1.23 


-0.28 


8.5E42 


0.19 


0.24 


10791. 


56200. 


1.90E13 


5.2E11 


-12.45 


-1.20 


-0.29 


9.2E41 


0.21 


0.22 


26411. 


27000. 


5.44E13 


1.8E11 


-13.21 


-1.02 


-0.24 


1.6E41 


0.23 


0.30 


44974. 


18700. 


8.92E13 


9.8E10 


-14.16 


-0.98 


-0.22 


1.2E41 


0.30 


0.36 


70104. 


14900. 


1.28E14 


7.5E10 


-14.36 


-0.95 


-0.14 


8.5E40 


0.33 


0.38 


99892. 


12100. 


1.70E14 


5.0E10 


-14.46 


-0.82 


-0.07 


8.3E40 


0.37 


0.43 



The explosion energy E^ n is not well determined by the absolute visual flux 
because of uncertainties in the distance and measurements, but can be estimated 
rather precisely by the rise time and the location of the first local minimum in 
V that depend on the shock travel time. Assuming that the time of the initial 
core collapse is given by the neutrino detections, a value of E^ n = 1.0 x 10 51 
erg can be ruled out. Since the shock travel time through the progenitor is 
essentially proportional to E®^, the energy must be larger than 1.20 x 10 51 erg. 
A value E^ in = 1.25 x 10 51 erg leads to a good agreement with the observations. 
Another, less stringent, upper limit can be derived from the upper limit set by 
Jones. SN1987A would have been seen by Jones if the shock travel time were 
about 30 minutes shorter. Therefore, E^i n must be < 1.6 x 10 51 erg. The hard 
UV fluxes are consistent with the observed high excitation lines of NV and OVI 
in the circumstellar ring (Fransson & Lundqvist, 1989). 

In the calculations just described no correction to the original temperature 
distribution given by flux-limited diffusion was made for the NLTE radiation 
field. This limitation was addressed by Mair et al. (1992) who used a gray 
radiation hydro code with adaptive mesh. For Arnett's 15 M model, Mair et 
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Figure 2. H a profiles with 
v s tat/v(Rph) [in %] of (solid), 
1 (long dashed- dotted), 5 
(dotted), 10 (short dashed), 20 
(short dashed dotted) and 50 
(long dashed) for a parameter- 
ized model with T e ff = 10, 000 
K, R ph = 1.3xl0 15 cm, v{R ph ) 
= 10,000 km s _1 and p oc r~ 7 
(Duschinger et al. 1995). 

al. found that E^in must be close to 1.3 x 10 51 erg and excluded 1.0 x 10 51 
erg. Ensman & Burrows (1992) and HE94 analyzed the explosion of Arnett's 
17 Mq model with an explosion energy of 1.0 x 10 51 erg. These analyses are 
similar to the approach just described using a combination of a gray radiation 
hydro code and NLTE-atmospheres (without time dependence). The explosion 
was followed up to day 4. As also stated in HE94, a problem seems to be that 
the photosphere recedes too rapidly and, consequently, the Doppler shifts of the 
absorption lines are too low after about day two. This suggests that a larger 
value of Eki n is required, as concluded above. 

4. The Photospheric Phase 

A variety of structures for the later photospheric phase were tested (e.g. H87, 
H88ab), based on explosions of models with a main sequence masses between 15 
to 20 M & provided by Weifi(private communication). The best agreement was 
found for 18 Mq. A statistical velocity field v stat was introduced to improve the 
agreement between observed and calculated spectra. The size of v sta t is small 
in comparison to v exp but larger than the sound speed, pa 15% of v exp during 
the first week and 5 % thereafter. This quantity should not be interpreted as 
microturbulence, but as a velocity field with a scale height smaller than the free 
mean path of the line photons (H87,H88a, see §6). With increasing v sta t, the 
features become broader, the emission component is more red shifted and the 
rise from the minima to the maxima of the P-Cygni profile is slower (Fig. 2). 

The luminosity was taken to be proportional to the observed bolometric 
LC. The choice of E^ n and the distance then determines the time evolution of 
the photospheric parameters such as R p h, T e ff, the radius Rhii of the outer 
boundary of the ionized hydrogen region, and the matter velocity at the photo- 
sphere (see Fig. 3). The spectral evolution and the colors are well reproduced 
as shown in Figures 4 and 5 (H87, H88ab). Spectral analysis of the photospheric 
phase gives E^ in = 1.3 x 10 51 erg, consistent with the shock break out and early 
LC analysis. The distance is discussed in §5. 

The difference between the spectral evolution of SN1987A and most other 
SN II can be understood as a consequence of the compact structure of the pro- 
genitor star which implies higher densities and steeper density gradients, and a 
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Figure 3. Distance R p h, 
T e ff, velocity v at R ph , 
the power law index h 
of the corresponding power 
law density slope and the 
distance Run up to which 
hydrogen is mainly ionized 
are given as a function of 
time. 
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low total luminosity For SN1987A, the recombination phase of the photosphere 
is encountered already after about 1 to 2 weeks. 

Early on, the location of Rph is strongly coupled to the expanding matter 
due to the high temperature and density at the photosphere. Consequently, the 
effective temperature T e ff drops very rapidly in the presence of a slowly varying 
total luminosity. This causes the rapid variations of the continuum flux in the 
UV, optical and IR during the first week. The smaller changes in the following 
months are due to the slow increase of R p h, due mainly to a geometrical dilution 
effect and to the recombination of H outside the radius Rhii- The decrease 
of Rph after maximum light is caused by the decreasing luminosity and by the 
higher helium abundances. Helium is mainly neutral and does not contribute 
to the electron density at the photosphere. Thus, the increase of the color 
temperature after maximum light is mainly a result of the increased helium 
abundance. The IR-excess can be understood as owing to free-free radiation 
and extension effects 10 to 20 %). No additional dust component is needed 
during the first few months at least in the nearer IR (H88ab). 

The observed frequency shifts of the absorption components of the Balmer 
lines are formed at different depths and hence probe the v/p structure (H91b). 
The good agreement between observed and calculated line shifts confirms the 
density profiles p(r). Lucy (1987) found similar density gradients from his analy- 
ses which are mainly based on the UV spectra. HE94 show fits of similar quality 
during the first 3 days with more self-consistent models (i.e. they calculate the 
luminosity from the radiation- hydro calculation). Their assumed p(r) is consis- 
tent with H87, H88, H91ab. At t=2d: p(R ph ) oc r~ 9 vs. oc r~ 10 in our work. 
Both in H87 and HE94 the p gradient changes over the line forming region by 
~ 20 %. Eastman & Kirshner (1989) used R p h, T^ejf an d p(r) oc r~ n as free 
parameters at each time. Although unrealistic, they used constant power law 
density profiles. They find gradients consistent with those given above (n(t=2 
d) = 9; n(t=10 d) = 5 to 7 vs. 10 and 7.2, respectively in H87). 

After one to two weeks, the Thomson scattering opacity declines due to 
recombination and heavy elements appear because weak lines are no longer 
"washed out." This spectral change does not imply a change in the composi- 
tion. Abundances less than solar are needed in order to reduce line blanketing 
by the metal lines. The abundances must be about one third of solar for nearly 
all elements (H88ab). This value closely resembles that derived for the LMC by 
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Figure 4. Spectra as 
observed by Menzies et al. 
(1987) February 25.9 and 
March 23rd, 1987, in com- 
parison with the synthetic 
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Dufour (1984). The exception in SN 1987A is the s-process elements Sc, Ti, V, 
Cr, Sr, Ba and Na. An overabundance relative to solar of Sc, V, Cr, Sr and Ti 
is needed to explain the slope of the spectrum below 5300 A. The strong line 
at about 6100 A can be attributed to Ba II (Williams 1988, H88ab) if Ba is 
increased by about 8-10 relative to LMC abundances. Other Ball features are 
present in the observed and calculated spectra at about 4900 A and 5000 A. 
They are consistent with the same overabundance; however blanketing by other 
lines prohibits a detailed analysis. Using Lucy's Monte Carlo code, Mazzali &: 
Chugai (1995) reanalyzed the spectra and found an overabundance for Ba of 
about 4 which is consistent within the uncertainties with other analyses. 

The evolving spectra are well reproduced by the atmosphere models during 
the first 5 months; however, later spectra (i.e. June, H88a) show insufficient 
blue shifts indicating too small a photospheric velocity. A power law density 
profile was introduced that yielded higher central densities than the original 
dynamical model. This is a crude way to allow for the entropy redistribution 
that must accompany mixing. At this epoch, excitation by 7-rays must also be 
taken into account because it strongly influences the temperature and ionization 
structure. The size of the 7-ray flux is assumed to be that given by the SMM 
satellite (Gehrels et al. 1989). 7-ray heating and inefficient cooling leads to high 
temperatures in the outer H-envelope and substantial emission by the cooling 
lines, e.g. Ca, despite their low abundance (Schwarz 1991, Li & McCray 1996). 

The calculated and observed spectrum on October 2 (Fig. 5) are in good 
agreement, if we choose roughly the same chemical abundances, but enrich He 
by a factor of 5. This implies that by this phase He-rich layers are seen. From 
the spectral analysis, the mass of the H-rich envelope is about 10 M (H88ab). 
A decrease of the H abundance is needed to explain the strong decrease of the 
continuum flux relative to the lines. The increase in the He abundance leads to 
a strong decrease of R p /t, since He is mainly neutral. Furthermore, this explains 
in a consistent way the increase of the color temperature given by B-V. At this 
time, hydrogen is still seen at the photosphere which shows an expansion velocity 
of the order of only 800 km s _1 . At the same time, CO bands in the IR indicate 
that the matter is strongly enriched with C and O up to radii that correspond 
to an expansion velocity of the order of 1800 km s _1 . These results strongly 
indicate that about 3M & of H have been mixed into very deep layers (H88ab; Li 
et al. 1993) of the ejecta and a moderate mass of C and O-rich matter, ~ 10~ 2 



Figure 5. Synthetic optical and f 
IR spectra corresponding to Octo- 
ber 2 nd , 1987, in comparison with 
observations in the optical (Oct. 
2 nd , Danziger, priv. comm.; left) 
and in the IR (ESO observations u 
between 1. and 6. Oct. 1987/ 
shifted by two units, lower right). 
C and O was assumed to be en- 
hanced to 20 * solar up to 1800 5 
km/sec. In addition, the theoreti- 
cal CO spectra are given for var- 
ious enhancement factors (lower 
right). Mixing out CO up to 3000 
km/sec would produce featureless 
CO-bands (Sharp & H6flich 1989). 
The feature at 2.38 fim cannot be 
attributed to CO+. 

to 1O _1 M0 have been mixed far out into the H-rich matter. This mixing can be 
understood by Rayleigh- Taylor instabilities during the explosion. Mixing was 
also indicated by the time dependence of the observed hard 7-rays, the detection 
of broad IR lines of Co and Fe, and the LCs ( e.g. Woosley, 1988). 

5. Distance to SN1987A 

Distances to many SN II have been derived by using the Baade-Wesselink 
method (e.g. SN 1969L: 12 ±4 Mpc and SN 1970G: 7±2 Mpc Kirshner & Kwan, 
1974; SN1979C: 23 ± 3 Mpc , Branch et al., 1981). In this method the ratios of 
the observed and predicted fluxes and spectral slopes at different wavelengths 
together with the velocities at different times are used to derive quantities such 
as the photospheric radius. In the past, these determinations were based on the 
assumption that the radiation of the supernova can be represented by a black 
body; however, the dilution of the radiation field and line blanketing are critical 
(Hoflich et al. 1986, Hershkowitz et al. 1986, H91c, Schmidt et al. 1994, Baron 
et al. 1995, Eastman et al. 1996). For high temperatures, the dilution factors 
are rather insensitive to the model parameters (10 to 25 %), but the absolute size 
given in the literature ranges from ss 0.3 (Eastman et al. 1996) to 0.5 (Baron et 
al. 1995) with H91c being somewhere in between. There is general agreement 
that the dilution factor rises steeply with declining T e ff during the recombina- 
tion phase of hydrogen (H91d, Schmidt et al. 1994). Typically, this phase is 
entered a few weeks after the explosion and is accompanied by the appearance 
of metal lines in the optical spectra. The critical T e ff for the onset of this rise 
depends sensitively on several parameters such as the structure and the He/H 
ratio. Unfortunately, uncertainties in the time since the explosion, At, enter 
the distance determination as At/t. This, together with the large sensitivity of 
the dilution factor as a function of T in the relevant range ~ 5000 K, poses a 
significant limitation on the accuracy of the Baade-Wesselink method; however 
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Figure 6. H a profile 

for different inclinations i 
for an oblate ellipsoid calcu- 
lated by a 3-D NLTE code. 
An axis ratio of 0.85 is as- 
sumed. Parameters repre- 
senting the atmosphere of 
SN1987A at about day 3 
after the explosion (T e ff= 
9000 K, R ph ,r sc =i = 8. x 10 14 
cm, p oc r~ 9 ). For details, 
see Hdflich et al. (1995) & 
Wang et al. (1997). 

high accuracy can be achieved if a detailed analysis of a broad range of data 
(including the time of the explosion) can be done. SN1987A can be regarded as 
the benchmark to demonstrate the "ideal" case. Based on the detailed spectral 
analysis for the first 200 days, (H88ab) derived a distance of 48 ± 4 kpc, well 
within the error bars of other groups (e.g. Chilukuri &: Wagoner 1988: 43 ± 4 
kpc; Schmutz et al. 1990: 46 kpc; Eastman & Kirshner 1989:49 ± 5 kpc). 

6. Deviations from Sphericity 

For SN1987A, polarization P of up to 0.5 % was observed early on. The level of 
P decreased rapidly over the next few weeks (e.g. Mendez et al. 1988). It was 
concluded that the envelope of SN1987A showed asphericities of about 10 to 20 
% and that Thomson scattering was the mechanism to produce P (Jeffery, 1991, 
Mendez et al. 1988, Hdflich et al. 1989, H91b). Whether P is caused by rapid 
rotation of the progenitor (Steinmetz & Hdflich 1992), by an asphericity of the 
SN Il-explosion mechanism (Yamada & Sato 1991), by an asymmetric luminos- 
ity input (H95), or dust in the surroundings (Wang & Wheeler 1996) is still a 
question under debate. There is growing evidence that, whatever the combina- 
tion of mechanisms, polarization is a rather common phenomenon among SNe II 
(Wang et al. 1995). The implications for our understanding of the explosion 
mechanism and progenitor evolution are obvious. There is also a direct impact 
on the use of SNe II as distance indicators because aspherical configurations also 
result in aspherical luminosities. An appreciable systematic bias will enter in 
statistical analyses when many SNe II are used to determine Hq because the 
directional redistribution of photons is different from the probability distribu- 
tion of the inclination angle seen by an observer (H91c). For a more detailed 
discussion of P see Wang, Wheeler & Hdflich (this volume). 

The finite polarization of SN II suggests a possible interpretation of the 
high statistical velocity deduced for SN 1987A as discussed in §4. An aspherical 
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configuration has some qualitatively similar effects on the spectra as a statistical 
velocity because both increase the spread in velocities observed at a given epoch, 
as can be seen by comparing Figures 2 and 6. Both cause the rise from the 
minima to the maxima to be less steep and the emission is shifted compared to 
a spherical configuration without v sta t (not shown in Fig. 6). In spite of the 
similarity, there does not exist a direct correspondence between the parameters 
for a spherical configuration v s t a t an d an aspherical configuration. More studies 
based on 3-D NLTE atmospheres are required (see Wang et al. 1997 for SNe la). 

7. Conclusions and Open Questions 

The peculiarities of the spectral and color evolution of SN1987A can be under- 
stood by the fact that the progenitor was a blue supergiant with its compact 
structure and low luminosity. Both from the early LC and the spectra we can 
conclude that E^ n is about 1.3 x 10 51 ± 0.3 erg based on stellar models with 
R p h = 3.10 12 cm. The hydrogen-rich envelope is about 9 — 11M . The uncer- 
tainty of the total mass of H is relatively small because of the weak dependence 
of the density at the photosphere on the atmosphere model parameters. Strong 
mixing of 3 Mq of H with the inner layers (down to 800 km s _1 ) and mixing of 
C/O rich layers out to about 1800 km s _1 was found. This can be understood by 
Rayleigh- Taylor instabilities during the explosion, but a full radiative transfer 
of a self-consistently mixed model has not been done. The metal abundances 
are consistent with the chemical composition of the LMC, but the s-process ele- 
ments are highly enriched. The abundance ratios of the s-process elements are in 
agreement with the s-process (Hashimoto et al. 1990), but the origin of the en- 
hancement remains unexplained. Polarization suggests asphericity in the outer 
layers of SN1987A of about 10 %. This may explain the high statistical velocity 
field needed to fine tune spectral fits during early epochs (H87); however, more 
detailed analysis is needed to test this suggestion. 

This research is supported by NSF Grant AST9528110. 
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